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Kaiser redshift-space distortion formula describes well the clustering of galaxies in redshift surveys on small 
scales, but there are numerous additional terms that arise on large scales. Some of these tenns can be described 
using Newtonian dynamics and have been discussed in the literature, while the others require proper general 
relativistic description that was only recently developed. Accounting for these terms in galaxy clustering is the 
first step toward tests of general relativity on horizon scales. The effects can be classified as two terms that 
represent the velocity and the gravitational potential contributions. Their amplitude is determined by effects 
such as the volume and luminosity distance fluctuation effects and the time evolution of galaxy number density 
and Hubble parameter. We compare the Newtonian approximation often used in the redshift-space distortion 
literature to the fully general relativistic equation, and show that Newtonian approximation accounts for most 
of the terms contributing to velocity effect. We perform a Fisher matrix analysis of detectability of these terms 
and show that in a single tracer survey they are completely undetectable. To detect these terms one must resort 
to the recently developed methods to reduce sampling variance and shot noise. We show that in an all-sky 
galaxy redshift survey at low redshift the velocity term can be measured at a few sigma if one can utilize halos 
of mass M > 10^^ Mq (this can increase to \0-a or more in some more optimistic scenarios), while 
the gravitational potential term itself can only be marginally detected. We also demonstrate that the general 
relativistic effect is not degenerate with the primordial non-Gaussian signature in galaxy bias, and the ability to 
detect primordial non-Gaussianity is little compromised. 
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I. INTRODUCTION 

In the past few decades galaxy redshift surveys have been 
one of the indispensable tools in cosmology, covering a pro- 
gressively larger fraction of the sky with increasing redshift 
depth. With the upcoming dark energy surveys this trend will 
continue in the future. However, despite the advance in obser- 
vational frontiers, there remained a few unanswered questions 
in the theoretical description of galaxy clustering. One is the 
issue of validity of the Kaiser formula, where density pertur- 
bation in redshift space is density perturbation in real space 
multiplied with a term that depends on the angle between the 
line-of-sight direction and the direction of the Fourier mode 
(we give a more detailed definition below) [1]. It has been 
well known (e.g., |01) that the simplest version omits some 
of the terms coming from the Jacobian of the transformation 
from real space to redshift space, terms of order v/Hr, where 
V is the velocity and r is the distance to the galaxy and H is the 
conformal Hubble parameter It is argued that these terms are 
potentially important, especially for large angular separations 
(e.g., see, ^3-;7::]), but most of the analyses so far have focused 
on effects in coiTelation function, without proper signal-to- 
noise analysis (see, however, [8J). 

A second, related issue, is whether the terms originally de- 
rived in the Newtonian approximation, get modified when 
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a proper general relativistic description is employed. On 
horizon scales, the standard Newtonian description naturally 
breaks down, and a choice of hypersurface of simultaneity be- 
comes an inevitable issue, demanding a fully relativistic treat- 
ment of galaxy clustering beyond the current Newtonian de- 
scription. In recent work [Io|, it is shown that a proper 
relativistic description can be easily obtained by following 
the observational procedure in constructing the galaxy fluc- 
tuation field and its statistics: We need to model observable 
quantities, rather than theoretically convenient but unobserv- 
able quantities, usually adopted in the standard method. While 
both the relativistic and the standard Newtonian descriptions 
are virtually indistinguishable in the Newtonian limit, they 
are substantially different on horizon scales, rendering galaxy 
clustering measurements a potential probe of general relativ- 
ity. 

The relativistic description of galaxy clustering includes 
two new terms that scale as velocity and gravitational poten- 
tial. Compared to the dominant density contribution, they are 
suppressed by H/k and (H/k)'^ and become important only 
on large scales, where the comoving wavevector amplitude 
is k. Consequently, the identification of these terms just by 
looking at the galaxy power spectrum is hampered because of 
sampling variance, and the general relativistic effects unac- 
counted in the standard Newtonian description may result in 
systematic errors less than l-cr for most of the volume avail- 
able at z < 3 in the standard power spectrum analysis [l^l- 
We revisit this calculation here, using a more realistic descrip- 
tion of these terms, but the basic conclusion remains the same. 
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In light of the fact that these corrections to the Kaiser 
formula provide a potential probe of the consistency of our 
model, including generic tests of general relativity on large 
scales, it is worth exploring if these terms can be observed. A 
new multi-tracer method 111 111 takes advantage of the fact that 
differently biased galaxies trace the same underlying matter 
distribution, and it can be used to cancel the randomness of 
the matter distribution in a single realization of the Universe, 
eliminating the sampling variance limitation. This method has 
been used in lfl2ll to investigate the velocity effects of fl fioll . 
noting that for any given Fourier mode the imaginary part of 
velocity couples to the real part of density and vice versa. 
Even with this novel technique, the expected detection level 
is low iflll . 

If sampling variance is eliminated, the dominant remaining 
source of error is shot noise, caused by the discrete nature of 
galaxies. Recently, a shot noise cancelling technique has been 
proposed 113 ] and investigated for detecting primordial non- 
Gaussianity UM in combination with the sampling variance 
cancelling technique. The basis of the method is that by using 
halo mass dependent weights one can approximate a halo field 
as the dark matter field and reduce the stochasticity between 
them. While this works best when comparing halos to dark 
matter, some shot noise cance lling ca n also be achieved by 
comparing halos to each other iImIIisII . This opens up a new 
opportunity to probe horizon scales and extract cosmological 
information with higher signal-to-noise ratio. 

In this work our primary goal is to explore the detectability 
of these effects in galaxy clustering on cosmological horizon 
scales using the galaxy power spectrum measurements with 
both single tracer technique and with the multi-tracer and shot 
noise cancelling techniques. In addition we clarify the rela- 
tion of the galaxy fluctuation field often used in the redshift- 
space distortion literature using Newtonian approximation to 
the fully general relativistic equation. Finally, we also inves- 
tigate the impact of the general relativistic effects in detecting 
the primordial non-Gaussian signature in galaxy bias. 

The organization of the paper is as follows. In Sec. |II] 
we present the full general relativistic description of galaxy 
clustering. In Sec. [nil we discuss the Newtonian correspon- 
dence and its relation to the redshift-space distortion literature. 
The multi-tracer and shot noise cancelling techniques are pre- 
sented in Sec.|IVl In Sec.[V] we provide the measurement sig- 
nificance of the general relativistic effects in the galaxy power 
spectrum. In Sec.|Vl] we extend our formalism to the primor- 
dial non-Gaussianity. Finally, we discuss the implication of 
our results in Sec. IVIII 



n. 



GENERAL RELATIVISTIC DESCRIPTION OF 
GALAXY CLUSTERING 



A full general relativistic descri ptio n of galaxy clustering 
is developed in [loll (see also Il6l - [l8ll ). Previously, we 
have adopted the simplest linear bias ansatz, in which the 
galaxy number density is just a function of the matter den- 



sity Ug = F[pm], both at the same spacetime.' However, this 
ansatz turns out to be rather restrictive, since the time evo- 
lution of the galaxy sample is entirely driven by the evolu- 
tion of the matter density oc {1 + z)^. Here we make one 
modification in the adopted linear bias ansatz by relaxing this 
assumption and providing more freedom, while keeping the 
locality. We allow the galaxy number density at the same mat- 
ter density to differ depending on its local history, as a local 
observer at the galaxy formation site is affected by its local 
matter density and the proper time to the linear order in per- 



turbation, i.e.. 



with tp being the proper time 



measured in the galaxy rest frame. Physically, the presence of 
long wavelength modes affects the local dynamics of galaxy 
formation by changing the local curvature and thus the ex- 
pansion rate, and these are modulated by the Laplacian of the 
comoving curvature and the proper time lil9l . Therefore, in 
addition to the contribution of the matter density fluctuation 
rasz = (5 — 3 (5z at the observed redshift z il[l3l,the physical 
number density of galaxies has additional contribution from 
the distortion between the observed redshift and the proper 
time, when expressed at the observed redshift: 



ng{z) [1 + & ((5 — 3 5z) — bt Sz^] 

ng{z) [1 + b 5y - e SZy] , 



(1) 



where the matter density fluctuation is S, the lapse 6z in the 
observed redshift z is defined as 1 + z = (1 + z)(l + Sz), and 
the homogeneous redshift parameter is related to the cosmic 
expansion factor as 1/a = 1 + z. The subscript v indicates 
quantities are evaluated in the dark matter comoving gauge. ^ 
The galaxy bias factor is 5 = d \n hg/d In pm\tp and the 
time evolution of the galaxy number density due to its local 
history is bt = dlnng/dln{l + z)\p^^. Therefore, the total 
time evolution of the mean galaxy number density is propor- 
tional to the evolution bias 1121 11. 



dlnn„ 



rfln(l 



3b + bt. 



(2) 



For galaxy samples with constant comoving number density, 
the evolution bias factor is e = 3 due to our use of the physical 
number density in Eq. (j2]i. This biasing scheme in Eq. ([T]i 
is consistent with 11614191 122tl . and our previous bias ansatz 
corresponds to bt = 0. 



' On large scales, a more general stochastic relation between the galaxy num- 
ber density and the matter density also reduces to the local form we adopted 
here [11]. As opposed to some confusion in literature, this biasing scheme 
is independent of whether galaxies are observed. 

^ Gauge-dependence arises, only when perturbation quantities are consid- 
ered. For example, the physical number density of galaxies is a well- 
defined scalar field without gauge ambiguity. However, when we split it 
into the homogeneous pait and the perturbation part, the correspondence of 
the physical quantity to the homogeneous part depends on the coordinate 
choice, and consequently the perturbation part becomes gauge-dependent 
(20[l . The comoving gauge is a choice of gauge conditions, in which v = 
or the off-diagonal component of the energy-momentum tensor is zero. No 
simple choice of gauge conditions corresponds to the observer's choice of 
coordinates {z, 0, 0). 
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Therefore, with this more physically motivated bias ansatz, 
the general relativistic description of the observed galaxy fluc- 
tuation is 



H- 



n 



-bp5VL-2K,. (3) 



Here the luminosity function slope of the source galaxy pop- 
ulation at the threshold is p. Note that this is the slope of the 
luminosity function, expressed in terms of absolute magnitude 
M — constant — 2.5 log L, hence the factor 5p instead of the 
factor 2a L, where is the slope of the luminosity function 
expressed in terms of luminosity L [see Eq. (|23]|1. In addition, 
the comoving line-of-sight distance to the observed redshift 
is r, the dimensionless fluctuation in the luminosity distance 
is SVl, the temporal and spatial metric perturbations are 
and tp^, the line-of-sight velocity is V, and the gauge-invariant 
radial displacement and lensing convergence are 6r^ and 
The subscript x indicates quantities are evaluated in the con- 
formal Newtonian gauge (also known as the zero-shear x = 
gauge). ^ We have ignored the vector and tensor contributions 
to (5obs in Eq. ©. 

We emphasize that compared to [l^ it is only the terms 
in the square bracket that are affected by the choice of linear 
bias ansatz, and Eq. ^ is consistent with lfl6l4l9ll22ll . Vari- 
ous other terms in Eq. (|3) arise due to the mismatch between 
the observed and the physical quantities. The radial and an- 
gular distortions are represented by 5r^ and fZ, and the dis- 
tortion in the observed redshift corresponds to the derivative 
term. The conversion of physical volume to comoving volume 
gives rise to the distortion 3 5z^, and the rest of the potential 
and velocity terms defines the local Lorentz frame, where the 
source galaxies are defined. For galaxy samples selected by 
the observed flux, additional contribution 5p SVl arises. Fur- 
ther discussion regarding this contribution is given later in this 
section. 

Since Eq. Q is written in terms of gauge-invariant vari- 
ables, it can be evaluated with any choice of gauge condi- 
tions. However, in evaluating Eq. (O, it proves convenient to 
use different combinations of gauge-invariant variables, rather 
than to choose one specific gauge condition and convert all the 
gauge-invariant variables in Eq. (|3]l to quantities in the chosen 
gauge. 



^ The displacement of the source galaxy position from the observed galaxy 
position is split into the radial and angular components, and these two parts 
are expressed in terms of gauge-invariant quantities, i.e., their forms differ 
in each gauge choice but give the same value regardless of gauge choice. 
These quantities are most conveniently expressed in the conformal Newto- 
nian gauge. In fioll the radial distortion is denoted as 57?.. Here we use a 
slightly different notation for the radial distortion 5r^ to avoid confusion 
with the dimensionless coefficient TZ. in Eq. )2U . 

The covariant derivative of the observer velocity is often decomposed into 
expansion, shear, rotation, and acceleration vectors |23|| . The shear compo- 
nent is proportional to x = + 7') with the metric convention in llOll . 
The conformal Newtonian gauge or the zero-shear gauge corresponds to 
the frame, in which there is no shear seen by a observer moving orthogonal 
to the constant-time hypersurface. 



For a presureless medium in a flat universe, the Einstein 
equations are (e.g., SEHH) 



fx 



•p'^ = Hay , 

and the conservation equations are 

v'^ + TLv^ = ka 
K = -3^; - kv. 



X 1 



X 1 



(4) 

(5) 
(6) 



(7) 
(8) 



where the prime is a derivative with respect to the confor- 
mal time T and the equations are in Fourier space. From 
these equations, it is well-known that the dark-matter comov- 
ing gauge {v = 0) is coincident with synchronous gauge 
[a^ = a„ + [av^y /ak, hence a„ = from Eg. 1©] and 
the comoving curvature is conserved ((/?', = 0) ll25ll26ll . Fur- 
thermore, since these gauge-invariant variables correspond to 
the usual Newtonian quantities iIzsIIztIi . hereafter we adopt a 
simple notation 6,, = dm, Vx = fx ~ ~'^x = to 

emphasize their connection. 

To facilitate evaluation of Eq. Q, we express the remaining 
gauge-invariant quantities in terms of 6m, v, and (p as 



dr J (27r)3 k 
5z = V + (j)+ [ df2(j>' 



Szy = Sz^ + 



Sr - 



d^k Hv 
I (2^ ~Y 

df2(t>, 







IC 



df 







(9) 
(10) 

(11) 
(12) 

(13) 



5Vl = (j) + V -^4r + I df-2(f>'- - I df2(f> 



dr 



Hr 



A</) + 0"-2-^ 
r V or 



St 

H± fe^ =-V- ^6' _ i±f ^ 
dz\H J H dr 

1 + zdH , 



(14) 



(15) 



where integration by part is performed in Eq. (Il4l and we have 
ignored quantities at origin that can be absorbed to the ob- 
served mean IsllT^l. The total derivative in Eq. (fT5t is with 
respect to the observed redshift z and it is related to the null 
geodesic as 



(16) 
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FIG. 1 : Redshift dependence of two dimensionless parameters 7?. and 
V in Eq. ( I21t . Non-vanisiiing values of 71 and V represent tiie gen- 
eral relativistic effects in galaxy clustering, each of which describes 
the contributions of the gravitational potential and the velocity to the 
observed galaxy fluctuation field. Three different curves represent 
galaxy samples in a volume-limited survey (p is constant) with three 
different limits Lt in luminosity threshold: a sample with low thresh- 
old Lt <^ L), (p — 0; dashed), a sample with no magnification bias 
Lt ~ L* (p = 0.4; solid), a sample at high luminosity tail Lt 3> i* 
(p — 1.5; dot-dashed). The evolution bias factor e = 3 is fixed in all 
cases, representing homogeneous galaxy samples (constant comov- 
ing number density) often constructed in large-scale galaxy surveys. 



The derivative term also appeared as a partial derivative in 
while we wrote it here as a total derivative to imply 
Eq. ( fTSI l. However, in a sense it is a partial derivative with re- 
spect to the observed redshift with other observable quantities 
{6, (/)) kept fixed. 

The observed galaxy fluctuation in Eq. ^ is the sum of the 
matter density dm, the gravitational potential (p, the line-of- 
sight velocity V, and other distortions such as 6z^, Sr^, and 
ST>L, and they are also a linear combination of 5„i, (f), and V 
with various prefactors and integrals in Eqs. (|9]— (fTSl), which 
in turn can be expressed in terms of the matter density dm- 
Using the Einstein equations, we have 



V 



3Hq r2„i 

k " ~ k 



(17) 
(18) 
(19) 



where the logarithmic growth rate is / = dliK^m/rflna, 
d/dr = Hd/dlna, and is the cosine angle between the 
line-of-sight direction and the wavevector. 

Before we compute the observed galaxy power spectrum, 
we make a further simplification by ignoring the projected 
quantities such as the gravitational lensing and the integrated 



Sachs-Wolfe contributions in Eqs. (fTOli— (fT4]i. which are im- 
portant only for the pure transverse modes (fc" =0) lfToll28ll . 
With this simplification and by using Eqs. (|9]|— (fTsll. the ob- 
served galaxy fluctuation in Eq. (O can be wiitten in Fourier 
space as 



d^k 
(2^ 



JNcwt 



VSm 

[k/nf 



USrn 

k/n 



(20) 

where the two dimensionless coefficients TZ and V are defined 

as 



V = ef--nm{z) 



e + .f- 



l + zdH 
H dz 



-(5p-2)(2-- 



(21) 



7^ = / 



l + zdH , / 1 

e ^ + (5p - 2) 1 - — 



H dz 



Ur 



and the standard Kaiser formula of the observed galaxy fluc- 
tuation is 



JNewt 



2 /, 



(22) 



which can be contrasted with (Jobs in Eq. (O. 

Apparent from their spatial dependence in Eq. (|20t . the co- 
efficients TZ and V originate from the velocity and gravita- 
tional potential. While Eq. ^ can be derived with the mini- 
mal assumption that the spacetime is described by a perturbed 
FLRW metric and photons follow geodesies, the coefficient V 
in Eq. (l2Tl i is obtained by applying the Einstein equations 
[Eqs. dUi— ©J. Given the observed e and p, the value and the 
functional form of V is, therefore, heavily dependent upon the 
general theory of relativity. 

While V is purely relativistic, some contributions to TZ 
may be considered non-relativistic, since they could be writ- 
ten down in Newtonian dynamics, simply as a coupling of 
velocity from the Doppler effect with the time evolution of 
the galaxy number density. However, different gravity mod- 
els will yield different values of TZ via the logarithmic growth 
rate / (see, e.g., ll29ll . and further discussion on this issue is 
presented in Sec.lIIIb. Therefore, measuring TZ and V is equiv- 
alent to a direct measurement of the relativistic contributions, 
and we collectively refer to the coefficients TZ and T^ as the 
general relativistic effects in the galaxy power spectrum. 

Figure [T] illustrates the redshift dependence of two coeffi- 
cients TZ and T^, in which the evolution bias factor is fixed 
e = 3 (constant comoving number density). The contribu- 
tions to these two coefficients arise from the volume and the 
source distortions that also involve the change in the observed 
redshift and the observed flux. Since the volume distortion in- 
volves 2Sr^/r in Eq. Q, both coefficients diverge at z — !> 
(r 0), unless the radial distortion of the volume effect is 
cancelled by the source effect (p = 0.4), leaving only the dis- 
tortion terms from the observed redshift (as we show below 
this cancellation is generic). With e = 3 and p = 0.4, the 
coefficient P in Eq. i2l[ nearly vanishes at low redshift. 

However, as we quantify in Sec. |V] even if the divergent 
term does not cancel out, since the survey volume decreases 
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faster, the diverging terms have neghgible impacts on the mea- 
surement significance. Furthermore, Figure [1] appears differ- 
ent from those obtained in ifisll . since they adopted the halo 
model to relate the evolution bias factor e to the galaxy bias 
factor b. Large-scale galaxy surveys show that these two pa- 
rameters for galaxy samples are independent of each other, 
and their relation is different from the halo model prediction 
(e.g., [30:]). 

Before we close this section, we discuss the fluctuation in 
the luminosity distance. Additional contribution 5p (52?^ to 
the observed galaxy fluctuation arises, because galaxy sam- 
ples are selected given an observational threshold in flux or 
equivalently a luminosity threshold Lt converted by using the 
observed redshift z. For flux-limited samples, the luminosity 
threshold is changing as a function of redshift, as only bright- 
est galaxies at high redshift can have observed fluxes large 
enough to be above the threshold in flux. Volume limited sam- 
ples are constructed by imposing a constant luminosity thresh- 
old Lt, up to a maximum redshift. Furthermore, the evolution 
bias factor e in this case is defined with respect to the number 
density of galaxies with L > Lt, in addition to other criteria 
such as color cuts.^ 

Fluctuations in luminosity distance bring galaxies around 
the threshold Lt into (or out of) the galaxy samples, and this 
effect provides the contribution 5p SV^ in Eq. OJ. The lu- 
minosity function slope of source galaxy populations at the 
threshold is defined as 



dM 



d log L 



where the absolute magnitude is related to the luminosity as 
M - Mi, = -2.5 log(i/L*) with pivot points Ah and L,^. 

The galaxy luminosity function is often described by the 
Schechter function llsill with a power-law slope as and an ex- 
ponential cut-off as 





L ' 


(£) '"P 





dL 



(24) 



and the number density given a threshold in luminosity is then 

ng{> Lt)^(^.T{as + l,Lt/L,) , (25) 

where T{a, x) is the incomplete Gamma function. Therefore, 
the luminosity function slope at threshold is 



p = OA 



Cts + l 



exp 



T{as + l,Lt/L^) 



(26) 



which can take values from zero to infinity, depending on the 
choice of Lt. 



A typical case of interest is a flux limited survey, which is 
dominated by L^-galaxies at the peak of the luminosity func- 
tion. For these galaxy samples p ~ 0.4 (solid), which can 
be obtained with Lt — and as = (or their variants). 
These galaxy samples have no diverging terms in TZ and V, 
because the volume distortion is balanced by the fluctuation in 
the luminosity distance. We may also assume volume-limited 
galaxy samples with constant p and consider two additional 
representative cases in Figure [T] First, galaxy samples at high 
luminosity tail (p = 1.5; dot-dashed), here taken as Lt = 3L^ 
and for which the Schechter function slope is as — — l.l(e.g., 
llsoll ). Last, we consider galaxy samples with sufficiently low 
threshold Lt <^ Li, {p ^ 0; dashed), which contain galaxies 
at low mass halos. 



III. NEWTONIAN CORRESPONDENCE AND 
REDSHIFT-SPACE DISTORTION 

In the cosmological context, full general relativistic equa- 
tions reduce to Newtonian equations on small scales, in which 
relativistic effects are negligible. This statement is born out 
by observations that matter density fluctuations 5 in various 
choices of gauge conditions become identical on small scales, 
except insofar as unusual gauge conditions are imposed such 
as the uniform density gauge, in which (5 = on all scales. 
However, as larger scale modes are considered, matter density 
fluctuations become increasingly different from one another, 
and the transition scale is largely set by horizon scales k = H. 

Even on these large scales, however, "Newtonian corre- 
spondence" exists in certain circumstances, in the sense that 
the matter density Sm, the velocity v, and the potential <f> fluc- 
tuations in Newtonian dynamics are identical to those quan- 
tities in fully relativistic dynamics with certain choices of 
gauge conditions ll25il27ll . In a flat universe with presureless 
medium, the Newtonian matter density 6m is identical to the 
comoving gauge matter density Sy, and the Newtonian veloc- 
ity V and potential cj) are identical to the conformal Newtonian 
gauge quantities and Lp^. Since this correspondence holds 
on all scales, numerical simulations properly capture large 
scale modes ll32l [33ll . apart from other technical issues such 
as finite box size. 

In redshift-space distortion literature iQl IH, a Newtonian 
calculation has been performed for obtaining the galaxy fluc- 
tuation field 5z in redshift-space that goes beyond the Kaiser 
formula in Eq. (l22l i. A different notational convention is of- 
ten adopted in the redshift-space distortion literature. Espe- 
cially, galaxy number densities Ug are expressed in comoving 
space. ^ The selection function a is defined with respect to the 



^ The luminosity function naturally evolves in time due to aging stars in 
galaxies, galaxy mergers, cosmic expansion, and so on. However, the con- 
dition that e = 3 and p is constant relies on the assumption that the shape 
of the luminosity function at constant Lt remains unchanged, while the 
physical number density ng(> Lt) is diluted as due to cosmic expansion. 



* Conversion of physical quantities to quantities in comoving space requires 
a redshift parameter. In observation, the observed redshift can be used 
without gauge ambiguity, but in general a redshift parameter is defined in 
conjunction with the expansion factor in a homogeneous universe. There- 
fore, it is a function of coordinate time and hence is gauge-dependent. Only 
in this section, did we use rig to refer to the galaxy number density in co- 
moving space, as there is no gauge ambiguity in the Newtonian limit. 
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comoving number density Ug of galaxies as 



din 7 



= 2 



rH 

l + z 



and the Une-of-sight velocity V is 



V 



l + z 
H 



V 



l + z 



(27) 



(28) 



where the last equality holds if we ignore potential contribu- 
tions to 5z^ in Eq. ( fTOl i. As a special case, the constant co- 
moving number density e = 3 corresponds to a = 2, and 
both of them are constant. 

The redshift-space distance s at the observed redshift z is 
then related to the "real" distance r as 



dz' 

IT 



l + z 
H 



r + V , 



(29) 



where the first equality is exact to the linear order in pertur- 
bation, while the second equality neglects potential contribu- 
tions to 6z^ From the conservation of total number of ob- 
served galaxies nz{s)d^s = nj.{r)d^r, one derives the rela- 
tion between the redshift-space and the real-space fluctuations 



1 + ^. 



ng{s) 



dh 



d^r 



s^ng{s) 



dV 

dr 



(30) 

and to the linear order in perturbations the galaxy fluctuation 
in redshift-space is then 



Sz ^ b S„,, - 



d_ 

dr 



V 



(31) 



l + zdV 2V l + zdH^^ 

This equation is known as the complete formula for the ob- 
served galaxy fluctuation in the redshift-space distortion lit- 
erature, while only the first two terms constituting the Kaiser 
formula |[l|] are often used. The additional terms come from 
the Jacobian of the transformation from real space to redshift 
space and also from the galaxy number density evaluated at 
the observed redshift. 

Even with the knowledge of the Newtonian correspon- 
dence, there are fewer terms in Eq. (|3ll compared to the full 
relativistic formula Sobs in Eq. (|3]l, and those terms account for 
missing physics in the derivation. Apparently ignored are the 
fluctuation SVl in the luminosity distance (effectively p = 0) 
and the lensing contribution that gives rise to distortions be- 
tween the observed and the source angular positions. Note 
that this by itself is an important effect: in a generic flux 
limited survey with p = 0.4 the potentially divergent term 



^ Since the real distance r takes a gauge-dependent redshift parameter z as 
an argument, we left the gauge choice of 5z unspecified, l + z = (1 + 
z){l + Sz). 



is exactly cancelled by the luminosity fluctuation effect, 
which is caused by the fact that a galaxy moved to a redshift 
closer to the observer will have its flux slightly smaller than 
what the redshift distance predicts, so it may not enter into a 
flux threshold, which compensates for the volume effect in the 
p = 0.4 case. 

However, once we account for luminosity threshold effects, 
the Newtonian calculation can fully reproduce the velocity 
terms in Eq. While the velocity terms receive a rela- 
tivistic contribution lfl2ll from the gradient of the potential in 
Eq. (fTsT i, it is cancelled by the time derivative of the veloc- 
ity via the conservation (Euler) equation in Eq. (|7). There- 
fore, the functional form of V is generic in Eq. ([3]l, same in 
all gravity theories (the conservation equations should hold in 
other theories of modified gravity, as it simply states that the 
energy-momentum is locally conserved ll34ll ). However, as we 
scale the velocity terms with the matter density Sm by using 
the Poisson equation (and also the conservation equation), the 
value of TZ itself will be different in other gravity theories than 
general relativity, and its measurements can test general rela- 
tivity, although this kind of tests can be performed by using 
the redshift-space distortion term in (5Ncwt- 

With respect to the relativistic contributions, adding the 
lensing contribution to the conservation relation in Eq. dSTJ is 
still not enough to recover the relativistic formula in Eq. (O. 
While the lensing contribution accounts for angular distor- 
tions, there exist a radial distortion in the source position, 
the Sachs-Wolfe effect Issll in the observed redshift, and fi- 
nally the difference in the observer and the galaxy rest frames. 
Equation ( |3TT i can be obtained from the relativistic formula in 
Eq. ^ by ignoring potential contributions, assuming p = 0, 
and identifying the matter density and the line-of-sight veloc- 
ity with those in the comoving and the conformal Newtonian 
gauges, respectively. However, the validity of the Newtonian 
description on large scales can only be judged retroactively, 
after the relativistic description is derived. We emphasize 
again that a fully relativistic treatment is required for estimat- 
ing P. 



IV. MULTI-TRACER SHOT NOISE CANCELLING 
TECHNIQUE 

We consider multiple galaxy samples with different bias 
factors for measuring the general relativistic effects in the 
galaxy power spectrum. Using a vector notation, the observed 
galaxy fluctuation fields can be written as 



If obs 



bo + ff^ll- 



cvV 

WW 



c-rR. 



k/n 



Sv 



= b{k, /ifc) +£ 



(32) 



where bo, I, £ are the hnearbias, the multi-dimensional iden- 
tity, and the residual-noise field vectors. By definition the 
noise-field is independent of the matter fluctuation (eSy) = 0, 
and the square bracket in Eq. (l32i defines the effective bias 
vector b. We will adopt a plane parallel approximation for the 
power spectrum analysis, meaning there is only one angle /ife 
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between the Fourier mode and the line-of-sight direction we 
need to consider The corrections to this approximation are 
expected to be small More generally, the effects consid- 
ered here are different from the plane parallel approximation 
and can be considered separately. 

The coefficients V and TZ in Eq. (ISTT i are also generalized 
to the multi-tracer case as 



-(5p 



7^ = / 



e 



- 2/) [2 - 
l + zdH 



1 



l + zdH 
H dz 



(33) 



(5p - 2/) 1 - 



We introduced two new parameters cn and c-p to generalize 
the measurement significance of the coefficients TZ and V to 
the case of multiple galaxy samples — they are ck = cp = 1 
in general relativity, and measurements of these two parame- 
ters amount to the measurement significance of the two vec- 
tors 'K and V. 

In order to assess our ability to measure the general rela- 
tivistic effects in the galaxy power spectrum, we employ the 
Fisher information matrix, and the likelihood of the measure- 
ments is 



1 



(27r)W/2 (detC)i/5 



cxp 



1 



-1 jfobs 
^GR 



(34) 



where the co variance matrix is C = (^gr ^qr ) = bb^ Pm + 
£, the shot noise matrix is 5 = {££^), and the matter power 
spectrum in the comoving gauge is Pm (fc). Since the observed 
galaxy fluctuation fields are constructed to have a vanishing 
mean (^q*^) = 0, the Fisher information matrix is 



- Tr [C-iQC-iCj] , (35) 



with two measurement significance parameters Oi — c-p, c-n 
and two nuisance parameter vectors 6i = e, p. The covari- 
ance matrix with subscript is Ci ~ dC/dOi. 

The inverse covariance matrix of the multi-tracer field and 
the derivative of the covariance matrix are 



£-^bWS-^P„ 



C, = 



1 



(36) 



= hb^ + bb] P 



where a = b^ S^^bP^ , bi = db/dOi , and we ignored the 
derivative of the shot noise matrix. With the inverse covari- 
ance matrix, we have 




10" 10'2 
minimum halo mass (h^'Mg) 



FIG. 2: Number density and average bias of halos above the mini- 
mum mass at different redshift slices. Since the multi-tracer method 
utilizes all the halos of mass above the minimum mass, a large num- 
ber of halos are required to achieve sufficiently low minimum mass. 



To further simplify the equation, we define two more coeffi- 
cients 



A - b^S-'b,P^ , 
7„ = Rc (bl5-^6,) P, 
and the various terms in Eq. 

b^C-'bP„, = 

Rc(blC-\)p,n - 



(39) 



b]£-% 



b\S-^b, 



Prr 



are 

a 
1 + a 

ft 



(40) 



1 



1^J 



13* + (3Wi 



2 (1 + a) 

The Fisher information matrix is therefore 

Re (ft/3, - aft/3;) 



F 



1 + a 



7ij 



(l + a)2 



(41) 



C-iQC-iCj- = P^i^C-^bblc-^bb] + 
+C-^b,b^C-^bb] + C-^b,b^C-^bjb'' 
and the Fisher information matrix is 

F,, = {b^C-'b)Re{blC-%)pl 

+Re [(fotc-ifo,) (btc-ib,)] p2 



^bblc-^bjb^ 



(37) 



This formula is a straightforward extension of the Fisher ma- 
trix in 111 41 with the effective bias vector b being a complex 
vector. The imaginary part arises solely from the 7?.-term in 
Eq. (|32] i and its derivative. 

For comparison, we also consider measurements by using a 
single galaxy sample. The formalism presented in Sec. llVlis 
valid for a single tracer, in which vector quantities reduce to 
scalar quantities. The Fisher information matrix for a single 
tracer is 



(38) 



2 Re(ft)Re(ft 
(l + a)2 



(42) 
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where we used the relation 



(43) 



only valid for a single tracer Note that in the single tracer we 
are not sensitive to the correlation between real and imaginary 
part of the mode, which has the dominant contribution to the 
signal-to-noise ratio of TZ in multi-tracer method [12]. This is 
only true in the plane parallel approximation. 



V. MEASUREMENT SIGNIFICANCE 

For definiteness, we consider full sky surveys with three 
different redshift ranges and adopt a set of cosmological pa- 
rameters consistent with the WMAP7 results Jsl!]. Given 
the survey volume V, the Fisher matrix is summed over 
the Fourier volume, where = 27r/y^/'^ and fci^ax = 

(we clarify the dependence of the measurement 
significance on our choice of fcmin and k^a.^). To model the 
Fisher matrix parameters a, jij, we adopt the halo model 
description in lfT4l[l5ll : It has been well tested against a suite 
of A'^-body simulations with Gaussian and non-Gaussian ini- 
tial conditions. 

We assume that the galaxy samples are constructed to have 
a constant comoving number density (e = 37) in a volume- 
limited survey (constant p = pi). While uncertainties in e 
and p can propagate to the measurement uncertainties in V 
and TZ, we focus on how well future surveys can measure V 
and TZ, assuming there are no uncertainties in theoretical pre- 
dictions of P and TZ. Given the current measurement uncer- 
tainties in e and p |30'], the uncertainties in theoretical pre- 
dictions are very small and will be smaller in future surveys. 
Figure [2] shows the number density and average bias of ha- 
los of mass above the minimum mass il/min at different red- 
shift slices. As we are interested in applying the multi-tracer 
method with sufficiently low minimum mass to enhance the 
measurement significance of the general relativistic effects, 
we consider two cases for the luminosity function slope at the 
threshold: p ~ (sufficiently low threshold Lt <C L^) and 
p = 0.4(Lt~L,). 



additional terms in Sz [Eq. dSTI )] used in the redshift-space 
distortion literature are all negligible in the standard power 
spectrum analysis (single tracer method). 

This result confirms the prediction in fu^] and extends the 
predictions to galaxy samples with different number density 
and bias. This conclusion is in apparent contradiction with 
i^, where correlation functions are shown with smaller errors 
than the size of the effects. However, these errors are obtained 
from simulations without taking into account the actual num- 
ber of modes in a realistic survey. Once this is taken into 
account there is no contradiction JH. The reason the wide- 
angle correction remains small is that there exists no velocity- 
density correlation ( ^ TZ) due to symmetry of pairs in Eq. (l20l i 
and the effect shows up only as a correction to the dominant 
contributions that are already accounted in the plane-parallel 
limit. 

The difference in the predicted measurement significance 
between the left {p = 0.0) and the right (p ~ 0.4) panels arise 
from the redshift dependence of the TZ and T^ values, as il- 
lustrated in Figure [T] The signals are computed at the mean 
redshift of each survey, and in the left panel {p = 0.0), the ab- 
solute values of TZ and T^ slowly decrease with redshift, while 
they remain nearly constant in the right panel (p = 0.4). Fur- 
thermore, even with higher suppression power of [k/Ti), it 
is generally easier to measure T' than TZ in the single tracer 
method for the case with p = 0.0 (left panel). Arising from 
the imaginary part in Eq. ( |32t , the 7^-term in the galaxy power 
spectrum is negligible, compared to the real part that includes 
the standard Newtonian term and the P-term, while the sen- 
sitivity to the 7-*-term comes from the cross-correlation of the 
standard Newtonian term and the P-term. In the left panel, 
the lack of sensitivity to the T'-term is due to the vanishing 
signal of T*. 

While the shot noise of massive halos is an obstacle for 
measuring the general relativistic effects, it is the cosmic vari- 
ance on large scales that fundamentally limits the measure- 
ment significance. Therefore, the measurement significance 
slowly increases as the minimum halo mass is lowered, but it 
quickly saturates given values of TZ and 7-". For measuring the 
general relativistic effects in the galaxy power spectrum, there 
is no further gain in constructing galaxy samples with large 
number density at a fixed survey volume. 



A. Single tracers 

First, we consider the prospect of measuring the general 
relativistic effects by using a single tracer Figure |3]shows the 
predicted measurement significance of TZ and for various 
survey redshift ranges. For all galaxy samples with different 
minimum mass (or different number density in Figure |2), the 
predicted measurement significance is very weak, indicating 
that substantial difficulty is present in measuring the general 
relativistic effects in the galaxy power spectrum for surveys 
at 2 < 3. This difficulty is simply due to the small number 
of large-scale modes that are sensitive to the general relativis- 
tic effects. Furthermore, the weak measurement significance 
of both TZ and T' means that compared to the standard New- 
tonian contribution (^Ncwt, the general relativistic effects or 



B. Multiple tracers 

This conclusion is, however, contingent upon the assump- 
tion that the power spectrum analysis is performed by using 
a single tracer, and the multi-tracer method can change the 
prospect of measuring the general relativistic effects in a dra- 
matic way. Figure |4] shows the predicted measurement sig- 
nificance of the general relativistic effects by taking full ad- 
vantage of the multi-tracer method. Compared to the single 
tracer case in Figure |3] there exist two key differences in the 
measurement significance derived by using the multi-tracer 
method. First, the measurement significance is greatly en- 
hanced in Figure|4]by eliminating the cosmic variance, which 
sets the fundamental limit in the single tracer method. Sec- 
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FIG. 3: Predicted measurement significance of general relativistic effects TZ (upper) and V (bottom) in tlie galaxy power spectrum obtained by 
using a single tracer. All halos of mass above minimum mass are lumped together to construct a single tracer. Four curves represent different 
survey redshift ranges with corresponding volume V — 2.5, 7.9, 59, 410 (/i~^Gpc)^. For the volume-limited sample (constant p) with 
constant comoving number density (e = 3), two galaxy samples are constructed to have p = (left) and p = 0.4 (right). No uncertainties in 
theoretical predictions are assumed. With the traditional power spectrum analysis (single tracer), it is difficult to measure the general relativistic 
effects at any meaningful significance. 



ond, a substantially larger measurement significance of the 7?,- 
term is obtained than that of the P-term by cross-correlating 
multiple galaxy samples and thereby isolating the imaginary 
term ||12!]. The method of measuring the imaginary part in 
the galaxy power spectrum of two tracers [12:] is fully im- 
plemented in our complex covariance matrix as off-diagonal 
terms and extended to the number of tracers larger than two. 
The result in lfl2ll would correspond to TZ/a ~ 1.8 at z < 1. 

In our most optimistic scenario, if we can utilize all halos of 
M > lO^^/i^^A/©, the velocity term TZ (solid) of the galaxy 
samples with p = (left panel) can be measured at more than 
10-CT confidence level at z < 1, while it is still difficult to 
detect the gravitational potential term V (solid). A significant 
detection of TZ can be made, even in surveys at low redshift 
(dotted and short-dashed), if halos of M < lO^^/i"^!/© can 
be used. At higher redshift z > 1, though the increase in the 
survey volume is partially cancelled by the lower abundance 
of halos at a fixed mass, a substantial improvement (dashed) 
for V can be achieved by going beyond z = 1, as the signal V 
increases with redshift (p = 0). 

However, the scenario above is not very realistic because 
of the p = assumption. In the right panels, we consider the 
galaxy samples with p = 0.4, of which the TZ and V values are 
nearly constant at all redshifts. The constant signals result in 



higher measurement significance for surveys with larger vol- 
ume at higher redshift. Compared to the case withp = 0.0, the 
measurement significance of TZ is smaller due to the smaller 
value of 7?. at z < 1, while that of T' is highly suppressed due 
to the vanishing value of T^. By using halos of mass slightly 
lower than Mq, a survey like the BOSS that covers 

redshift range z = 0.3 ~ 0.5 with a quarter of the sky can 
achieve a 1-ct detection of TZ, demonstrating the feasibility of 
the multi-tracer analysis in future surveys. 

A few caveats are in order First, we used the galaxy 
power spectrum in a flat-sky and counted the number of modes 
in computing the measurement significance. Calculations of 
wide-angle correlation functions ISfl shows that the effects of 
geometry are negligible, and a similar calculation can be per- 
formed for power spectrum measurements (in preparation). 
Furthermore, our calculation is rather insensitive to the mini- 
mum wavenumber fcmin- At sufficiently low Mmin, the multi- 
tracer method approaches the optimal case with dark matter 
density field, where the Fisher information can be approxi- 
mated as b'£~^b' Pm- In terms of spatial dependence alone 
in Eq. (l32T i. F-p (x J dk k'^{k^ /k"^) oc ln(fcniax/fcmin) and 
Ftc oc /dfc P{k"/P) oc (fc^iax ~ ^min) ^ith the spectral 
index ri ~ 1. Therefore, the dependence of the measurement 
significance on fc„iin is logarithmic for V and negligible for TZ, 
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FIG. 4: Predicted measurement significance of general relativistic effects TZ (upper) and V (bottom) in the galaxy power spectrum derived 
by using the multi-tracer method. All halos of mass above minimum mass are utilized to take advantage of the multi-tracer method |[TTl[l4ll . 
Various curves are in the same format as in Figure |3] Compared to Figure[3] the measurement significance is substantially enhanced by using 
the multi-tracer method. 



as fcmax ^ fcmin- In coiitrast, a significance enhancement in 
F-ji can be achieved by increasing the maximum wavenumber 
^max, although the gain is marginal for F-p. Finally, while 
any degeneracy with cosmological parameters in measuring 
the relativistic effects is largely eliminated due to the cancel- 
lation of the underlying matter distribution [HTlfTil . a proper 
quantification of the measurement significance requires con- 
siderations of uncertainties in theoretical predictions of P and 

7^. 



VI. PRIMORDIAL NON-GAUSSIANITY 

We extend our formalism to the primordial non-Gaussian 
signature in galaxy bias ifjTll and introduce additional param- 
eter /nl. Here we only consider the simplest local form of 
primordial non-Gaussianity to demonstrate how it can be im- 
plemented in the general relativistic description, and ignore 
scale-independent and scale-dependent corrections (see, e.g., 

MM)- 

The primordial non-Gaussian signature in galaxy bias can 
be readily implemented in our full general relativistic descrip- 
tion with the Gaussian bias factor in Eq. Q replaced by 



& + 3/NL(fc-l)^c 



where 5c is the linear overdensity for spherical collapse 
and is the transfer function of the curvature perturbation 
(see also ifisi [l9l I22II ). Equivalently, the primordial non- 
Gaussianity can be considered as additional relativistic con- 
tribution by replacing V in Eq. ( 1211 1 with 



(45) 



Figure |5] shows the predicted constraints on the primordial 
non-Gaussianity derived by accounting for the full general rel- 
ativistic effects. In obtaining the constraint (7/^^ on primor- 
dial non-Gaussianity, we set c-jz — c-p = 1 and marginalize 



(44) 



over e and p with priors (7,. = 0.1 and ap — 0.05 OOll . We em- 
phasize that e and p can be more accurately measured in obser- 
vations, further reducing their uncertainties. With the current 
uncertainties in e and p, the constraint (7^^^ (solid in Fig.|5]) is 
nearly identical to the unmarginalized constraint. The dashed 
curve shows that even with no priors on e and p, a /^^l is not 
inflated except in the regime with ct/^l 2, because TZ and 
V are affected simultaneously by e and p but only V by /nl. 
Furthermore, the unique dependence of /nl on 6 — 1 and 
in Eq. ( l44l i provides the multi-tracer method with more 
leverage to separate it from the general relativistic effect. It is 
the primordial gravitational potential at initial epoch, not the 
evolved gravitational potential at the observed redshift that the 
scale-dependent galaxy bias responds to, and the difference is 
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no prior, no slope 

z^l no prior with slope 

0.1 I — — — — — 

109 IQIO IQll IQIZ 1013 1014 

minimum halo mass (h^^M^) 

FIG. 5: Predicted constraints on the primordial non-Gaussianity /nl 
from galaxy power spectrum measurements. To facilitate the com- 
parison, the constraints on /nl are obtained by using the same sur- 
vey specifications as in IT4II : V ~ 50 (/i~^Gpc)^ centered at z = 1. 
We assume e = 3 and p — 0.4, and various curves show crf^^ with 
different priors on e and p (a^ — 0.1, ap — 0.05 jsoll for the solid 
curve). 

the transfer function T^p{k,z). At low redshift, the transfer 
function decays from unity on scales smaller than the horizon 
around the dark energy domination epoch, and there exists a 
factor two difference in r^(fc, z) at our adopted fcmax- 

Finally, we allow e and p to vary as a function of mass with 
two logarithmic slope parameters a,, and a^: 

e = eoI + aeln(M/Mo) , (46) 
p = po/ + apln(M/Afo) , 

with eo = 3, po = 0.4, ae = ap = 0, and Mq = 
lO^^/i^^A/0. The effects of and Up (dotted in Fig.|5]l are 
sufficiently different from that of /nl, and (T/„l asymptoti- 
cally reaches the floor set by the uncertainties in e and p. This 
demonstrates that the general relativistic effects in the galaxy 
power spectrum are not degenerate with the primordial non- 
Gaussian signature. However, if /nl were to be constrained 
below unity, similar precision needs to be achieved in predict- 
ing TZ and V. 

This requirement can be relaxed by increasing the max- 
imum wavenumber /cmax to exploit the unique dependence 
of /nl on T^p{k,z). With larger maximum wavenumber 
^max = O.l/iMpc"^, the overall uncertainties on /nl are re- 
duced by about 30% for all three cases in Figure |5] Another 
way is to construct galaxy samples by using independent mass 
estimates instead of observed flux, as p drops out in Eq. ( |20l i, 
further separating e from /nl- 



VII. DISCUSSION 

In this paper we explore the contributions to the redshift 
survey beyond the Kaiser approximation in the context of re- 
cently developed general relativistic analysis. We compare 
the results of this formalism to the previous analyses |j6l and 



show that these analyses ignore several terms such as lumi- 
nosity distance fluctuation and overestimate the significance 
of the effect. In addition, correlation function analyses pre- 
viously adopted are not optimal to assess the signal-to-noise 
ratio of these effects, in contrast to our power spectrum analy- 
sis. We find that these corrections beyond the Kaiser formula 
are not observable in a generic redshift survey using a single 
tracer, meaning that these effects do not have to be considered 
in a generic redshift survey. A caveat to this conclusion is that 
in this paper we have performed the analysis using the plane 
parallel approximation and do not consider large angle effects 
IE 01 J but we expect that these effects are equally small 
Their detectability will be addressed in a separate publication. 

Using the multi-tracer shot noise cancelling method the de- 
tection significance is increased, providing a unique oppor- 
tunity to test these effects, and general relativity in general, 
on horizon scales. Still, for realistic cases the detection sig- 
nificance is at the few sigma confidence level (in some more 
optimistic cases the detectability rises to 10 sigma or more), 
so this test of general relativity is of interest only when the de- 
viations from general relativity are significant on large scales. 
We have also shown how the primordial non-Gaussian effect 
in galaxy bias can be implemented in the full general rela- 
tivistic description, and we have argued that the ability to de- 
tect primordial non-Gaussianity is little compromised by the 
presence of general relativistic effects. 

Considering the fact that we perform a very large-scale 
analysis, our method of measuring the general relativistic ef- 
fects in galaxy clustering need not be restricted to spectro- 
scopic surveys. The use of photometric redshift measure- 
ments may not affect our results if the photo-z error is suf- 
ficiently small, e.g. an error of dz/(l + 2) ^ 0.03 corre- 
sponds to fc > 0.06/iMpc~^ at 2 = 1, which is larger than 
^max = 0.03/iMpc^^ we adopted here. This allows one in 
principle to extend the observed halo mass ranges to lower 
masses using a photometric survey and to take full advantage 
of the multi-tracer method. 

While we treated multiple galaxy samples as halos in mul- 
tiple mass bins and our method requires sufficiently low 
^J^min, the SDSS already measures galaxies well below AI = 
lO^^h^^ Mq (e.g., SDSS LI sample), and there exists numer- 
ous methods to remove satellite galaxies and isolate central 
galaxies. Furthermore, it is shown iffill that one needs a fairly 
large scatter in the mass-observable relation, a\nM = 0.5, 
to degrade the shot-noise suppression, and the scatter is less 
important at the low mass end. The reason for this insensi- 
tivity to the scatter is that in terms of weighting, the opti- 
mal weighting method puts more weight on massive halos, 
w{M) = M + Mq, where A/q is a constant and approximated 
as 3 times the minimum mass. We note however that our pre- 
diction is based on the halo model description of the shot noise 
matrix, which is tested only for halos at M > 10^^ /i^^ A/0, 
and our prediction at AI < 10^'^ AIq is an extrapolation. 

The bottom line of this paper is that the corrections beyond 
the Kaiser formula in redshift surveys are generally small and 
only detectable using very specialized techniques adopted in 
this paper. This is a good news for those analyzing generic 
redshift surveys since they do not have to consider them. Nev- 
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ertheless, the potential detectability of these terms gives rise 
to the prospect of testing general relativity in a regime previ- 
ously untested. Thus despite the small detection significance 
it is worth exploring these tests further to see if they can be of 
use in separating general relativity from some of its alterna- 
tives. 
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